Anomalous decays of 77' and rj into four pions 
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We calculate the branching ratios ol the yet unmeasured rj' decays into four pions, based on a 
combination ol chiral perturbation theory and vector-meson dominance. The decays rj' — > 2(7r+7r~) 
and ri' — > 7r^7r~27r° are P-wave dominated and can largely be thought to proceed via two p reso- 
nances; we predict branching fractions ol (1.0 ±0.3) x 10~* and (2.4 ±0.7) x 10~*, respectively, not 
much lower than the current experimental upper limits. The decays rj' ^tv^ and 77 — ^ 47r°, in con- 
trast, are D-wave driven as long as conservation ol CP symmetry is assumed, and are significantly 
further suppressed; any experimental evidence for the decay rj — >■ 4^° could almost certainly be 
interpreted as a signal of CP violation. We also calculate the CP-violating amplitudes for ri' — >■ 47r'' 
and r) — >■ 47r'' induced by the QCD ^-term. 
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I. INTRODUCTION 

Processes in low-energy QCD that involve an odd num- 
ber of (pseudo-)Goldstone bosons (and possibly pho- 
tons), which are, therefore, of odd intrinsic parity, are 
thought to be governed by the Wess-Zumino-Witten 
(WZW) term |lj via chiral anomalies. While the so- 
called triangle anomaly is well tested in processes such 
as tt", 77 — >■ 77, and the box anomaly contributes e.g. to 
77r — ;> TTTT and 77 — >■ ttttj, the pentagon anomaly remains 
more elusive; the simplest possible process that is usu- 
ally cited is K'^K~ — >■ 7r"'"7r~7r'', which however has not 
been experimentally tested yet, and is likely to be subject 
to large corrections to the chiral-limit amplitude that is 
dictated by the WZW term. 

A different set of processes involving five light pseu- 
doscalars is the four-pion decays of rj and rj' . Experimen- 
tal information about these is scarce: only upper limits 
on branching ratios exist 0; however, this may change 
in the near future for at least some of the possible final 
states with the advent of high-statistics rj' experiments 
such as BES-III, WASA-at-COSY, ELSA, CB-at-MAMI- 
C, CLAS at Jefferson Lab, etc. We are only aware of 
one previous theoretical calculation of these decays, per- 
formed in the framework of a quark model whose 
partial width predictions, however, have in the mean- 
time been ruled out by the experimental upper limits, at 
least for the channel ?/ 2(7r+7r^). 

In principle, the decays rj' — 47r, in contradistinction 
to many other 77' decay channels, seem not terribly forbid- 
den by approximate symmetries: they are neither isospin 



forbidden, nor required to proceed via electromagnetic 
interactions. The reaction rj — > 47r, in contrast, is es- 
sentially suppressed by tiny phase space: only the decay 
into 47r° is kinematically allowed {M^-4:Mt,o = 7.9 MeV, 
Mr,-2{M^±+M^o) = -1.2MeV). Furthermore, the fact 
that anomalous amplitudes always involve the totally an- 
tisymmetric tensor ef^i,ap can be used to show that no two 
pseudoscalars are allowed to be in a relative S wave: as- 
suming they were, this would reduce the five-point func- 
tion PPPPP effectively to a four-point function SPPP 
(where S stands for a scalar and P for a pseudoscalar), in 
which there are no four independent vectors left to con- 
tract the e tensor with. The decays 77' — >■ 2(7r"'"7r~) and 
77' — > 7r+7r~27r° can therefore be expected to be P-wave 
dominated. As furthermore Bose symmetry forbids two 
neutral pions to be in an odd partial wave, 77' — >■ Att'^ 
and 77 — )■ 47r° even require all tt'^ to be at least in rel- 
ative D waves [1]. This, combined with the tiny phase 
space available, leads to the notion of 7/ — 47r" being CP 
forbidden [2, 5, 6], although strictly speaking it is only 
S-wave CP forbidden. 

The outline of the article is as follows. We begin by dis- 
cussing the two decay channels with charged pions in the 
final state, 77' -J> 2(7r+7r~) and 77' 7r+7r^27r°, in SecHH 
There, we calculate the corresponding decay amplitudes 
at leading nonvanishing order in the chiral expansion, 
saturate the appearing low-energy constants by vector- 
meson contributions, and calculate the corresponding 
branching ratios. In Sec. IIIIl we then construct a CP- 
conserving (D-wave) decay mechanism for 7/, r/ — )■ Att'^ 
and determine the resulting branching fractions, before 
discussing the CP-violating (S-wave) 77, 77' — >■ 47r° decay 
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as induced by the QCD 6'-term in Sec. IIVI Finally, we 
summarize and conclude. The Appendices contain tech- 
nical details on four-particle phase space integration as 
well as on a (suppressed) tensor-meson mechanism for 
77, 77' —J- 47r". 
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II. T]' 2(7r+7r-) AND t]' 7r+7r-27r° 

A. Chiral perturbation theory 

We wish to calculate the leading (nontrivial) chiral con- 
tribution to the anomalous decays 

T]' 7r+(pi)7r"(p2)7r+(p3)7r"(p4) , 

77' ^ 7r+(pi)^°(p2)7r-(p3)vr*'(p4) ■ (1) 




FIG. 1: Feynman diagrams contributing to 77' 2(7r'*'7r ) 
(and similarly to 77' 7r+7r~27r°) at 0{p^). The thick dot in 
the right diagram denotes a vertex from -Cq^^. 



The amplitudes can be written in terms of the invari- 
ant variables Sij = (pi + pj)'^, j — 1, ... ,4, which are that only contribute for the singlet field 770: 
subject to the constraint 

S12 + Sl3 + Sl4 + S23 + S24 + S34 = M^, + Mil (2) 



(in the isospin limit of equal pion masses). The five- 
meson vertices of the WZW term can be deduced from 
the Lagrangian 
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where is the number of colors and will be taken to be 3 
in this paper, Ft^ — 92.2 MeV is the pion decay constant, 
and (. . .) denotes the trace in flavor space. For simplicity, 
we refrain from spelling out the WZW term in its full, 
chirally invariant form. Furthermore, 
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We assume a simple, one-angle 7777' mixing scheme, 

I77) = cos6'p|778) - sin6'p|77o) , 
It/') =sin6'p|778) +cos6lp|77o) , 



(4) 



(5) 



and use the standard mixing angle 9p — arcsin(— 1/3) « 
— 19.5°. As we are going to present what in some sense 
corresponds to a leading-order calculation of the decay 
amplitudes, we regard the more elaborate two- angle mix- 
ing schemes Q as beyond the scope of this study; we 
expect the error made thereby to be covered by our gen- 
erous final uncertainty estimate. 

The flavor structure of Eq. is such that there are 
no direct contributions to 77, 77' — > 47r, and the decay am- 
plitudes vanish at leading order (in the anomalous sec- 
tor) 0{p^). Nonvanishing contributions only occur at 
O (p^ ) , where the amplitudes are given by sums of (kaon) 
loops and counterterm contributions from the 0{p^) La- 
grangian of odd intrinsic parity @, see Fig. [T] Only 
two different structures (esc CY , C^) remain when ex- 
ternal currents are switched off. Ref. Q only considers 
the Goldstone boson octet; we add terms oc C}^, 



-^e^,a.0{h^''[u^,u''u'']){u^) + ... , (6) 

with the usual chiral vielbein = i{u^d^u — udf^u^) 
(neglecting external currents), u = exp(7(p/2_F'7r), /i^jy — 
V^u^ + \7^Uf, with V^X = df,X + [r^,X] and = 
^{u^df^u + ud^v)) (neglecting again external currents). 
Furthermore, we use %_ = u^xm^ — wx^u, where x — 
25 diag(mu, TTT-d, 771^) -t- . . . contains the quark mass ma- 
trix and B is related to the quark condensate according 
to B = ~{qq)/F^. The decay amplitudes at 0{p^) take 
the compact forms 



AVO/8 TI'+TT TT+TT ) = -^(770/8 ^ TT+TtV Tt") 
■P1P2P3P4 [-^0/8(512) + -7^0/8(534) 
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Ms) = -16x/2 (c^/ifi) - Ci^'^(m) J s , 
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1, 
(7) 



with the scale-dependent renormalized low-energy con- 
stants C^"^ (/i) and CI2 (m) ■ There are no loop contribu- 
tions to the 770 amplitudes at this order, since at 0{p'^) 
the anomalous five-pseudoscalar term ([3]) (the left ver- 
tex of the loop diagram in Fig. [ij contributes only to 
the octet case. Equation ([7]) is scale-independent with 
the P function for C}^''(^) obtained in Ref. Q, if we 
demand C12 to have the same infinite part and scale de- 
pendence as Ci2- A numerical estimate for the finite 
part Cjf''(Mp) will be obtained by resonance saturation 
through vector-meson contributions. 
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B. Resonance saturation from hidden local 
symmetry 

Resonance saturation for the 0{p^) chiral Lagrangian 
of odd intrinsic parity has been studied in great gener- 
ahty recently in Ref. Here however we opt for the 
simpler, but on the other hand more predictive hidden- 
local-symmetry scheme [lOl - fTsj , which has the additional 
advantage of having been tested phenomenologically in 
great detail [l4| . 

In the framework of hidden local symmetry (HLS), 
there are four additional terms involving vector-meson 
fields, with coefficients Ci (i = 1, . . . , 4)jjn addition to the 
WZ W action for anomalous processes Ell ; as already 
noted in Ref. [151], only three independent combinations 
of these contribute to low-energy amplitudes at 0[p^). 
HLS amplitudes for any given anomalous process con- 
tain two kinds of contributions: contact terms and reso- 
nance exchange terms. The contact terms have the same 
form as those derived from the WZW action, but with a 
modified coefficient (see below); the gauge-invariant con- 
struction of the HLS Lagrangian density guarantees that 
the additional, Ci-dependent contributions will be can- 
celed by vector-meson exchange in the low-energy limit. 
In the following, we again for simplicity reasons refrain 
from properly defining all the HLS Lagrangian terms in 
their chirally invariant forms, but only quote the terms 
relevant for vertices of five pseudoscalars; the full La- 
grangians can be retrieved e.g. from Refs. [12, 13]. 

The contact terms for five-pseudoscalar vertices can be 
read off from the Lagrangian 



^HLS _ Nce^^a^ 
2407r2F5 



1 - — (C1-C2) 



(8) 

The low-energy limit of the vector-meson-exchange con- 
tribution can be obtained by integrating out the heavy 
fields: substituting the leading-order equation of motion 
of the vector-meson fields 



1 



(9) 



where g is the universal vector-meson coupling constant, 
into the HLS Lagrangians [H, [T^ 



Hvvp 
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327r2i^3 



(10) 



where the vector-meson nonet (with ideal mixing) is de- 
fined as 
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we find 



Ncjci - C2) 
1287r2F5 



p {ipd^'ipd''ipd"ipd^ip) , (12) 



which exactly cancels the second term inside the square 
brackets in Eq. ^ . 

If we extend the equation of motion Eq. ^ to next- 
to-leading order in the derivative expansion. 



1 



d' 



(13) 



where My is the vector-meson mass, we can derive the 
vector-meson contribution to the five-meson vertices at 
O(p^). Inserting Eq. ^ into Eq. we find 



C 



(6) _ Ncjci - C2 + C3) 



X {d^d''ip[dxip,d''>pd°'Lpdl^Lp\ 
- 2d'^ipd^ipd''ipd'^ipd^ip) . 



(14) 



The first term is exactly of the form of the Lagrangian 
term cx in Eq. ([6]). For the second term, we use 
the equation of motion for the Goldstone bosons, which, 
neglecting higher orders in the fields, reads (compare e.g. 
Ref. 0) 



(15) 



so we also identify a vector-meson contribution to 
and C]^^ . Our results read altogether 



1287r2M2 



2Ci^'-(Mv 



(16) 



where we have indicated the conventional assumption of 
the resonance-saturation hypothesis to be valid roughly 
at the resonance scale, /i = Aly (which in the following 
we will identify with the mass of the p, Mp = 775.5 MeV). 
The numerical values of the HLS couplin g co nstants are 
often taken to be given by ci — C2 ~ C3 w 1 [10] , fairly con- 
sistent with more elaborate phenomenological fits that 
yield ci - C2 = 1.21, C3 = 0.93 [11]. 

In principle, this completes the task to provide the nec- 
essary input for an evaluation of the chiral representation 
of the decay amplitude, Eq. ([T)). We observe, however, 
the following. First, evaluating the slope of (the largely 
linear function) J-b{s) in Eq. ([7]) with this input (using 
j;^^(0) = l/(967r2M2.)), we find 



87r2(47rF^)2 x ^"8(0) 



3(ci - C2 -I- C3 



(471^^2 

2Af2 



1.21og^). (17) 



Numerically, the first term is about 6.7 x (ci — C2 -l-C3)/2, 
and the second is 0.1. Hence, at the scale p. = Mp, the 
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slope is entirely dominated by the vector-meson contri- 
bution, and the kaon loops are negligible. 

Second, the maximal value for the kinematical invari- 
ants in 77' — >■ 47r allowed by phase space is y^SiJ < 
M^' — 2A/^ fa 680 MeV, therefore replacing the p propaga- 
tor by its leading linear approximation is not phenomeno- 
logically reliable. Even deviations induced by the finite 
width of the p, Tp = 149.1 MeV, will be clearly visible. In 
the following, we will therefore use the full vector-meson- 
exchange amplitudes as derived from the HLS formalism, 
with the p-meson propagators including the width, which 
in addition is expected to be a very good estimate of the 
higher-order pairwise P-wave interaction of the pions in 
the final state (of course neglecting any crossed-channel 
effects). They are given by 

■^viVS Tr+TT^TT+TT^) = —r=Av{r]0 Tr+TT^TT+TT^) 

v2 

= -Avivs^ 7r+7r"7r^7r") = ^Avivo 7r+7r°7r^7r°) 

v2 
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is the inverse p propagator, and we have neglected 
the width term in the transformation from Eq. ()18p to 
Eq. ([T^ in order to demonstrate the correct chiral di- 
mension O(p^) of the vector-meson contribution explic- 
itly. Expanding the resonance propagators in Eq. ([19]) 
and comparing to Eq. ([7]) easily leads back to the cou- 
pling constant estimate for C|^'' found on the Lagrangian 
level in Eq. p^ . 

At this point, we can try to answer the introductory 
question on which parts of the WZW anomaly action — 
triangle, box, or pentagon — the decays rj' — )■ 2(7r+7r^) 
and 77' — >■ 7r+7r^27r*' yield information. As the pentagon 
anomaly only enters via the kaon-loop contributions, we 
have found above that its significance for the decays un- 
der investigation here is negligible; the vector-meson con- 
tributions are derived from the triangle and box-anomaly 



terms, see Eq. ([10)) . As the phenomenological values of 
the HLS coupling constants suggest ci — C2 — C3 <C 2c3, 
the box anomaly yields the lesser part of the two, and 
the decays are dominated by the triangle-anomaly term. 



C. Branching ratios 



We calculate the partial widths of the decays rj' 
r~7r"'"7r~ and tj' — 7r"'"7r^7r~7r" using 

^ 4^) = J \A{rf ^ 4^)pd$4 



(21) 



where the evaluation of the four-particle phase space 
$4 is discussed in detail in Appendix |X] S* is a sym- 
metry factor — 5* = 4 for the 2(7r+7r~) final state, and 
S — 2 for the 7r+7r~27r° one. Note that with the relation 
A{r]o — >• 4tt) = y/2A{r]s — >■ 47r) and the standard mixing 
according to Eq. ([5|) , we have A{r]' — > An) = A{ris — )• An) . 
To obtain branching ratios, we normalize the partial 
widths by the total width of the rj' as quoted^ by the 
particle data group, T^/ = (0.199 ± 0.009) MeV 0]. Note 
that by using the most precise single measurement of 
this width alone, E^/ = (0.226 ± 0.017 ± 0.014) MeV [3, 
our predictions for the branching fractions would be re- 
duced by more than 10%. Given the observation of 
Eq. ([r7[). we neglect the kaon-loop contributions alto- 
gether and evaluate the matrix elements using Eq. ()18p . 
In order to account for trivial isospin-breaking effects due 
to phase space corrections, we calculate the branching 
ratio for rj' n^n'^n^n'^ using an average pion mass 
= (Mjr+ -f M^o)/2, while we employ the charged pion 
mass for the decay into four charged pions. All results 
are first quoted as a function of the coupling constants 
ci — C2 and C3, before inserting two sets of values: (i) 
ci - C2 = C3 = 1, and (u) ci - C2 = 1.21, C3 = 0.93 [li]. 
We refrain from employing the errors given in the fits 
in Ref. IJj: the uncertainties in the HLS coupling con- 
stants are well below what we estimate to be the overall 
uncertainty of our prediction. The results are 

B{t]' 2{n+n-)) 

= 0.15 (ci - ca)^ + 0.47 (ci - C2)c3 + 0.37 4 

= {1.0, 1.1} X 10"^ , 
6(77' ^ n+n-2n°) 

0.35 (ci - C2)^ + 1.09 (ci - C2)c3 + 0.87 

{2.3, 2.5} X 10"^ . 



X 10-^ 
(22) 



X 10" 



(23) 



We therefore find that the uncertainties due to the HLS 
coupling constants are small. We wish to point out that 
although nn P-wave dynamics are usually well approx- 
imated by the p resonance, and crossed-channel effects 
are expected to occur rather at the 10% level (as inferred 
from studies of decays such as w — )• 37r, (/)—)• Stt [13]), the 
present study in some sense still amounts to a leading- 
order calculation: SU(3)-breaking effects of the order of 
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Fjj/Ftt ~ 1.3 [III may occur, and in the treatment of 
the rj' (?7o), we have imphcitly evoked the l/N^ expan- 
sion. We therefore deem a generic uncertainty of 30% 
reahstic, and quote our predictions accordingly as 



B{ri' ^ 2(7r+7r-)) = (1.0 ± 0.3) x 10"'' 
B(tj' 7r+7r"27r°) = (2.4 ± 0.7) x lO"'' 



(24) 



These are to be compared to the current experimental 
upper limits 0, [l^ 



> 2(7r+7r-)) < 2.4 x lO"'' 
Tr+7r-27r") < 2.6 x 10"^ 



(25) 



hence signals of these decays ought to be within reach of 
modern high-statistics experiments soon. 



A. Pion-loop contribution 



Our decay mechanism for 77, ry' 47r° is built on the 
observation that there is a specific diagrammatic contri- 
bution that we can easily calculate, and that, in partic- 
ular, comprises the complete leading contribution to the 
imaginary part of the decay amplitude. This is given by 
7r+7r~ intermediate states, and hence harks back to the 
results of the previous section. As argued above, it ap- 
pears at chiral 0{p^'^): rj^/g, — )■ 7r+7r~27r*' as calculated in 
Eq. ^ to 0{p^), followed by rescattering tt+tt^ — > 7r°7r°, 
where the S wave does not contribute, and D and higher 
partial waves start to appear at 0{p'^) [l^l; see Fig. [5] 
for illustration. We calculate this first in the follow- 
ing approximation: given the numerical dominance of 



III. r,, r,' A-K° 

As we have mentioned in the Introduction, the P-wave 
mechanism described in the previous section, proceed- 
ing essentially via two p intermediate resonances, cannot 
contribute to the 47r'' final states. In fact, we can show 
that the D-wave characteristic of 77, 77' — >■ An^ suppresses 
these decays to 0{p^^) in chiral power counting, that is 
to the level of three loops in the anomalous sector. This 
is, in particular, due to the flavor and isospin structure of 
the anomaly, which does not contain five-meson vertices 
including 27r° at leading order (©(p"')), and to the chiral 
structure of meson-meson scattering amplitudes, which 
only allows for S and P waves at tree level {0{p'^)). As 
a complete three-loop calculation would be a formidable 
task and is certainly beyond the scope of our exploratory 
study, we instead consider the decay mechanisms shown 
in Fig. [2] As shown in Appendix [b1 the contribution from 
two /2 mesons is negligible in comparison to the pion 
loop. We therefore focus on the pion loop as shown in 
the left panel of Fig. [2l It represents a decay mechanism 
that, we believe, ought to capture at least the correct 
order of magnitude of the corresponding partial width. 




4tv". The 




FIG. 2: Left: Pion-loop contribution to rj, 77' 
black circle denotes an effective local ri, rj — > 7r^7r~27r'' cou- 
pling at C'(p^), the black square an effective local D-wave tttt 
scattering vertex at 0(p*). Right: 77, rj — )■ 47r'' through two 
intermediate /2 mesons. 



the counterterm contribution in Eq. ([T]), the amplitudes 



J^o/8(^) ^-re taken to be linear, J-qi^{s) 



-^0/8(^)'*' 



glecting tiny curvature effects from the kaon loops; and 
we approximate tttt rescattering by a phenomenological 
D wave, thus improving on the leading chiral represen- 
tation, but neglecting G and higher partial waves. We 
find 



^(778 ^ 47r") = ■^-^('yo -> 47r") = ^ — ^-^-iip5'p^p^p^|5(si2, S23, S34; S13) 

+ ^/(Sl2, Si4, S23, S34; S24) - Q{siz, S23, Sl4, S24; S12) - ^?(si3, Si4, S23, S24; S34) 

- ^/(Sl2, S24, Sl3, S34; S14) - ^(Sl2, Sl3, S24, S34; S23)} , 

g(v,w,x,y-s) = 3(5-4Af2)2 |(^-4A4) •^--(«) 

- 2(s2 _ lOsAf 2 + 30Af4) fi + 1 log -i- f f! - \sMl + XhU^ 



167r2 ° p I 167r2 ^15 
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<J^\ 1- 



4Af2 



L = 



d-4 



1 



167r2 \ d-4: 



2^1 E- l-log47r) 



(26) 



t2(s) is the partial wave of angular momentum 1 = 2 for 
the appropriate isospin quantum number /; the expres- 
sion 164(s)(s - 4M2)-2 = a| + 0{s - AM^), with the 
D-wave scattering length 03, is therefore finite at thresh- 
old. Note furthermore that = C(p^) in chiral counting, 
such that the chiral order of Eq. (|26p is indeed ©(p^"). 
L contains the infinite part of the divergent loop dia- 
gram in the usual way, using dimensional regularization. 
Of course, this individual loop contribution is both di- 
vergent and scale-dependent: only the imaginary part is 
complete (to this order) and in that sense well-defined 
and finite. We display the full expression here as we will 
use the scale dependence as a rough independent consis- 
tency check below. 

Without the knowledge of counterterms of an order 
as high as ©(p^^), one cannot make a quantitative pre- 
diction using the loop amplitude derived in the above. 
Hence, we have to resort to a certain phenomenological 
representation. The imaginary part of Eq. (|26p. which is 
complete at 0(p^°) as mentioned, is used to establish a 
connection to a one-/2 exchange in the s channel. Note 
that the /2(1270) exchange dominates the available tttt 
scattering phase shifts in the I = 0, £ ^ 2 channel, see 
e.g. Ref. gll. 

We will proceed to estimate the full D-wave tttt rescat- 
tering contribution as follows. Neglecting again any 
crossed-channel effects, rescattering of two pious can be 
summed by the Omnes factor. 



exp 



Sj{z)dz 



(27) 



where 5l is the tttt scattering phase shifts in the channel 
with isospin / and angular momentum (.. Near threshold, 
its imaginary part can be approximated as 

luvniis) « 5l{s){l + 0{a^))^atl{s){l + 0{a^)] 

(neglecting the shift from unity in i7(4M2), which is jus- 
tified in the D wave for our intended accuracy), while in 
the approximation of a phase being dominated by a nar- 
row resonance of mass M and width F, the Omnes factor 
is given by 



M^exp {i5l{s)) 



AMI 

— —\ r 

-ami) 



M / s- 



(29) 



Despite D-wave scattering near threshold not being dom- 
inated by the /2( 1270),^ we still use Eq. ([28]) to invoke 
the /2. This is because at somewhat higher energies, the 
/ = TTTT D wave dominates over the 1 = 2 component 
and will be well approximated by the /2(1270) resonance. 
One may wonder whether, in particular, for 77 — > An'^, 
which stays close to tttt threshold throughout the allowed 
phase space, this approximation may not lead to sizeable 
errors. We have checked for the numerical results for the 
branching fraction discussed below that, employing the 
full Omnes function according to Eg^ ipTl) with the phase 
parameterization provided in Ref. [23], the branching ra- 
tio changes by about 10%, well below the accuracy we can 
aim for here. On the other hand, within the ry' — > An^ 
decay, we stay sufficiently far below the resonance en- 
ergy that the phase of the D wave can still be neglected. 
With the correspondence between Eqs. ((28)) and (|29|) . we 
conclude that the /2(1270) contribution to the amplitude 
can be estimated as 

1 



_ Ncjci - C2 + C3) e^^g/g 



IJ' V a P 
P1P2P3P4 



247r2 

1^^/2(312, S23, Sl4, S34; Sl3, S24) 

— Gf2 (Sl3, S23, Sl4, S24; S12, S34) 

- ^/2(si2,S24,Sl3,S34;Sl4,S23)} , 



Qf2iv,w,x,y;s,t) 



V ~ w — X + y 



Af2 



M2 

12 



Mj 



M^f 

J2 



Mj 

J2 



(30) 



neglecting for simplicity the width of the /2, which is 
justified in the kinematic regime accessible in rj' — Att'^ . 
Note that, due to the special symmetry of the ampli- 
tude, Eq. (j30p can be rewritten identically by employing 
a "twice-subtracted" version of the resonance term, i.e. 
replacing Gf^ G'J^ , 



GfAv,w,x,y]s,t) = 



M^M^ 

P 12 



Mj 

J2 



M? -t 

J2 



(31) 



^ It is dominated by the low-energy constant £2 from the 0{p'^) 
Lagrangian ,20], or by t-channel vector- meson exchange in the 
spirit of resonance saturation |22ll . 
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TT 



,0 



471" via in- 



FIG. 3: Pion-loop contribution to r], 77' 
termediate states; see the vector-meson dominated amplitude 
discussed in Sec. [TTl The black square denotes an effective 
local D-wave tttt scattering vertex at C'(p^). 



the finite counterterm by 



IGtt 



xAf4 



0.22 



(33) 



where we have used 0° 
-3^-4 " 



1.75 X 10-3M-4, a\ = 0.17 x 
10~'^M~"^ [2^. In other words, the scale dependence sug- 
gests the order of magnitude of our counterterm estimate 
using /2 saturation to be reasonable. 



which makes the correct chiral dimension of the reso- 
nance contribution manifest. 

As a rough final consistency check, we compare the 
order of magnitude of a chiral counterterm induced by the 
Ji exchange, see Eq. ([5T|) in the low-energy limit s, t 
Mj^ , with the scale running of such a counterterm as 
necessitated by the log fi dependence in Eq. If we 

only retain the scattering lengths in the D-wave partial 
waves, the relevant part to be compared to Eq. (PT|) (that 
does not cancel in the full amplitude) is 



-(02-02) 



87r2 



(32) 



Comparing the numerical prefactors, we find that the 
scale dependence is suppressed versus the estimate for 



B. Pion-loop contribution improved: including 
vector propagators 



We have seen in Sec. |TI] on the P-wave dominated, 
(partially) charged four-pion final states that the lead- 
ing approximation in an expansion of the p meson prop- 
agators is not a sufficient description of these decays, 
given the available phase space in 77' decays. With the 
Vo/s ^ 7r'^7r"27r*^ transitions entering the decay mecha- 
nism for rjQ /8 — >■ Att'^ as described in the previous section, 
this deficit would be fully inherited in our estimate of the 
all-neutral final states. In fact, the imaginary part of the 
corresponding diagram including the full p propagators, 
see Fig. [21 can even be calculated exactly, using Cutkosky 
rules; however, the resulting expressions are extremely 
involved and not very illuminating. It turns out, though, 
that the main effects of the not-so-large vector-meson 
mass can be approximated by the following expression 
for the imaginary part: 



^/2 Stt V3i^5 



PlP2?'3P4{(ci - C2 - C3) Init;f(si2,S23,Sl4,S34;Sl3) 



-I- Imtjf (si2, si4, S23, S34; S24) - Im^^f (si3, S23, Si4, S24; S12) - Im tjf (si3, S14, S23, S24; S34) 

- Imtjf (si2, S24, Si3, S34; S14) - Imtyf (si2, Si3, S24, S34; S23) 

2C3 Imt/^(si2, S23,Si4, S34; S13) + lmQ^{si2, Sl4, S23,S34; S24) " Ini^?2(^13, ^23, Sl4, S24; S12) 

- ImCJ^(si3, si4, S23, S24; S34) - Im^^(si2, S24, S13, S34; S14) - ImCJ^(si2, S13, S24, S34; S23) | 



(AP^^^iv + w))' (Af2_l(^ + y))2 



Im g^{v,w,x,y; s) 

(M"^ (v — W — X + y) — Vy + Wx) , . o , <7 ^ v\ 

Img^iv,w,x,y;s) = -^^^ JJ ' i4{s) - tl{s))- + O [a') 

{M^p-h{v + w)) {Mj-\{x + y)) 37r 



(34) 



We find, furthermore, that the neglected terms indicated as ©(cr^) are also suppressed in inverse powers of Mp, 
starting at 0{M~^) compared to the leading terms of 0{M~^) in the above. Numerically, the indicated higher-order 
corrections in are found to be small, less than about 10% all over phase space. However, the corrections by the 
remnants of the p propagators are large compared to the limit Adp — >■ cx), given the available phase space and the 
high power of these propagators in the denominator. Using the same trick as in the previous section to transform the 
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imaginary part into an estimate for the whole (resonance-dominated) partial wave via the Omnes function, we arrive 
at 



Aim 



1 



4^") = -^^(r;o -> 47r") 



c ^^va/S 



24^2 



PlP2P3Pi[i(^l -C2-C3) i(si2,S23,Si4,S34;Sl3) 



+ S/2^i(si2, Sl4, S23, S34; S24) - ^^/2^i(si3, S23, Si4, S24; S12) - Gf^^iisi3,Sii, S23, 824, 834) 

- Gf^^i{si2,S24, Sl3, S34; S14) - Gf^ iisi2, Sl3, ^24, ^34; S23) 

2C3 Gf^^2i^l2,S23: Sl4, S34; S13) + Gf^^2i^l2, Su, S23, S34; S24) " Gf^^2i^l3: S23, Sl4, S24; S12) 

- 2(^13' *14, S23, S24; S34) - ^/2_2(*12> «24, S13, S34; S14) - ^/^^ 3 ('5l2 i ^13 , ^24, ^34; S23) | , 

M^{v-w) Mj{x-y) 
M^(M2(u - w - a; + y) - + wx) 



£2 



/2 



(Af2 -\{v + w)Y (M2 - i(a; + y))^ ^^/^ 



(35) 



Note that this result is far from the one-/2 dominance 
estimate, with a /2 coupling constant oc Mr^ as the 
previous section suggested. Expanding Ec^. (1351) simul- 
taneously around the limits Mp — oo, Mf^ — oo, the 
leading term (corresponding to chiral dimension Oip^'^)) 
is not dominated by terms of 0{M~'^MJ^), but also con- 
tains other terms of 0{M-'^MJ^) and 0{M-^). In other 
words, Eq. pop is numerically no reasonable approxima- 
tion to Eq. (|35|) even for the decay 77 — >■ 47r°, with its tiny 
phase space available. 



C. Branching ratios 

We calculate the partial width using Eq. ([2T|) with the 
symmetry factor 5 = 4!. Note again that A{ri' 47r") — 
A(ri^ ~^ 47r°), assuming standard mixing. We employ 
the amplitude as given in Eq. as our "best guess" 
for an estimate of the branching fraction. With the same 
numerical input as in Sec. Ill CI (except using the neutral 
pion mass everywhere), we find 

6(77' -> 47r") 

= [0.4 (ci - C2f + 1.6 (ci - C2)c3 + 1.7 c2] X 10-8 
= {3.7, 3.9} X 10"* , (36) 

for the two sets of coupling constants c^. Note that the 
use of the amplitude ([50]) leads to a branching fraction of 
the order of 4 x 10"^^, i.e. almost 3 orders of magnitude 
smaller. 

We can trivially also calculate the branching fraction 
for rj — )• 47r° , the only 77 — 47r decays that is kinematically 
allowed. We again employ the amplitude and note 
that mixing according to Eq. ([5]) suggests A{ri — > 47r") — 
\/2A{iis, — >■ 47r'^). Normalized to the total width of the ry. 



r,, = (1.30 ± 0.07) keV [1], we find 
6(77 ^ Ati°) 

= [0.4 (ci - C2f + 1.1 (ci - C2)C3 1.0 c2] X 10-30 
= {2.4, 2.6} X 10-3" , (37) 

in other words, the D-wave characteristic of the decay 
combined with tiny phase space leads to an enormous 
suppression of the CP- allowed r\ — \t\P decay. We again 
compare these estimates to the available experimental 
upper limits d, , 

Scxp(r/' ^ 4^0) < 5 X 10-4 , 
Scxp(»7 ^ 47r") < 6.9 X lO"'^ ; (38) 

further improvements of these experimental upper limits 
are planned (see e.g. Ref. for r\ — )■ 47r0). In this case, 
our predictions are smaller than those by several orders 
of magnitude. 

The uncertainties of Eqs. and ([57)) are hard to 
assess. The generic SU(3) and 1/A^c error of about 30% 
assumed in Sec. Ill Cl is probably too small, as here, we do 
not even have a complete leading-order calculation at our 
disposal. We therefore rather assume these numbers to 
be the correct orders of magnitude, without quantifying 
the uncertainty of the prediction any further. 

IV. CP- VIOLATING r], t)' 4Tr° DECAYS 

Given the smallness of the branching fractions pre- 
dicted for 77' 47r0, 77 — >■ 47r° via a D-wave dominated, 
CP-conserving decay mechanism in the previous section, 
it is desirable to compare these numbers with possible 
CP- violating contributions that may, on the other hand, 
avoid the huge angular-momentum suppression. One 
such CP-violating mechanism that is expected to affect 
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27| 



strong-interaction processes is induced by the so-called 6- 
term, an additional term in the QCD Lagrangian necessi- 
tated for the solution of the U(l)yi problem. The 6'-term 
violates P and CP symmetry and may induce observ- 
able symmetry-violating effects, in particular, in flavor- 
conserving processes. Its effective-Lagrangian treatment 
includes a term that can be rewritten as (see Ref. 
and references therein) 

^« = *'°^V^-^^-i°nd^, 

U^u^^ exp (^) , (39) 

which, in addition to the well-known 77 — >■ 27r ampli- 
tude [13, H^, also induces a CP- violating 77 — )■ 47r ampli- 
tude, 

Acpivs 47r°) = -^Acpivo -> 47r°) 

= AcpW 4.°) = ^-^CPiv ^ 4.°) = . 

(40) 

We will use M^^ « M,y for numerical evaluation. The 
fact that this amplitude is a constant makes the phase 
space integration almost trivial, with the results for the 
branching fractions 



^4^0) = 5x10-^ 



I CPV 



47r°) = 9 X 10^2 X el 



(41) 



We remark that we do not consider the branching ratio 
estimate for 77' — > Att^ in Eq. (PT|) reliable in any sense: 
given the available phase space and the possibility of 
strong S-wave tttt final-state interactions, it could easily 
be enhanced by an order of magnitude. Were ^0 a quan- 
tity of natural size, Eq. (|^T|) would demonstrate the en- 
hancement of the C P- violating S-wave mechanism com- 
pared to the CP-conserving D-wave one, see Eqs. 
and (157)1 . With current limits on the QCD vacuum angle 
derived from neutron electric dipole moment measure- 
ments, ^0 ^ 10~^^ [2^, these branching fractions are 
already bound beyond anything measurable; however, 
we note that for 77 47r°, the suppression of the CP- 
conserving D-wave mechanism, see Eq. ([37]), is so strong 
that it is even smaller than the CP-violating (S-wave) 
one in Eq. (|41l) if the current bounds are inserted for 6q. 



V. SUMMARY AND CONCLUSIONS 

In this article, we have calculated the branching frac- 
tions of the rj and 77' decays into four pions. These pro- 
cesses of odd intrinsic parity are anomalous, and — as long 
as CP symmetry is assumed to be conserved — forbid the 
pions to be in relative S-waves. We organize the am- 
plitudes according to chiral power-counting rules, and 



find the leading contributions to the 77' decay amplitudes 
with charged pions in the final state at 0{p^). Utiliz- 
ing the framework of hidden local symmetry for vector 
mesons, we assume that vector- meson exchange saturates 
the 0{p^) low-energy constants, and find that the (P- 
wave) decay amplitude is entirely governed by p interme- 
diate states. The dominant contribution is hence given 
by the triangle anomaly via 77' — )■ pp (with numerically 
subleading box terms), not by the pentagon anomaly. In 
this way, the branching fractions for 77' — >■ 2(77+77") and 
77' — > 7r+7r" are predicted to be 

6(77' ^ 2(7r+7r")) = (1.0 ± 0.3) x 10"^ , 
B{ii 7r+7r"27r°) = (2.4 ± 0.7) x 10"^ , (42) 

respectively. The former is only a factor of 2 smaller 
than the current experimental upper limit, so should 
be testable in the near future with the modern high- 
statistics facilities. 

Predictions for the decays into four neutral pions are 
much more difficult, as Bose symmetry requires them 
to emerge in relative D-waves (assuming CP conserva- 
tion), suppressing the amplitudes to 0{p^'^) in chiral 
power counting. We here do not even obtain the full 
leading-order amplitudes, as these would require a three- 
loop calculation. We estimate the decay via a charged- 
pion-loop contribution with D-wave pion-pion charge- 
exchange rescattering; an alternative mechanism through 
two /2 mesons is found to be completely negligible in 
comparison, based on an estimate of the tensor-tensor- 
pseudoscalar coupling constant in the framework of QCD 
sum rules. Because of these phenomenological approxi- 
mations, the CP-conserving branching ratios thus ob- 
tained, 



6(77' ^ 47r°) - 4 X IQ-s ^ 
B(77 ^ 47r°) ~ 3 X IQ-^o ^ 



(43) 



should only be taken as order-of-magnitude estimates. 
It thus turns out that the CP-conserving decay width 
of 77 — > 47r° is so small that any signal to be observed 
would indicate CP- violating physics. For the latter, we 
calculate one specific example using the QCD 0-term. 
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Appendix A: Four-body phase space integration 

The n-body phase space is defined as 
d^n{P;Pi, ■ ■ ■ ,Pn) 

Using the recursive relation Q 

d$„(P;pi, . . . ^ d<i>j{q;pi,...,pj) 

X (i$„_j+i(P;9,Pj+i,...,p„)— , (A2) 

we have 

d$4(-P;Pi,---,P4) 

dq^ dk^ 

= d^2 iq;pi,P2)d^2ik;P3, ^4)^*2 (P; g, fc) — — 

ZTT ZTT 

/> A I— 7113—7114 rM — si2 

/ mi +7712 ^ m3+m4 



1 

37r2)4M 



X / dnidn'^dn\pi\\p'^\\q\ , 



(A3) 



where TOj , i = 1 ... 4 arc the masses associated with the 
final-state particles of momentum , M is the mass of the 
decaying particle, S12 = g^, S34 = fc^. dfJ^ = dc/jjdcos^i 
is the solid angle of particle 1 in the center-of-mass frame 
(cmf) of particles 1 and 2, dVl'^ is the solid angle of par- 
ticle 3 in the cmf of 3 and 4, and dVt is the solid angle of 
the 1, 2 system in the rest frame of the decaying particle. 
The three-momenta are given by 



\Pi\ 



|g| = 



2M 



IP3I 



Ai/2(s34,m2,m|) 



2x/S34 



(A4) 



with the usual Kallen function A(x, y, z) = + — 

2(xy + xz + yz). 

Denoting quantities in the cmf of 1 and 2 (3 and 4) by 
*('), one can relate them with those in the rest frame of 
the decay particle by Lorentz transformation. Explicitly, 



Pt - {1i2{pT 



P2 



{712 (p? 



,7l2(/3l2P?* 
,7l2(/3l2P2* 



-p*l\\),p*ll} , 
-P*l\\),-Pl±} 



Ps = {734(P§' +/334 •P3),734(/334P§' +P3i|).P3±} , 

{734(rf' - /334 • P^), 734(/334rf' - Psll), -PsA , 



P4 



(A5) 



where ^12 = q/g*^ (/334 = k/k^) is the velocity of the 
1,2 (3,4) system in the rest frame of the decay particle, 
and 712(34) (1 - /3?2(34))"^^^- Moreover, pt||(_L) are 
the components of pi parallel (perpendicular) to q, and 
P3||(±) components of pjj parallel (perpendicular) 




FIG. 4: The solid angles of particle 1 (3) in the cmf of 1 and 2 
(3 and 4) . &13 is the angle between the momentum of particle 
1 in the cmf of 1 and 2 and the momentum of particle 3 in 
the cmf of 3 and 4. 



to k. One can define 61 as the angle between the direc- 
tions of q and pi, and 6*3 as the one between k ^ q and 
P3. The angle between pi and p'^, 613, is related to the 
solid angles il^ and ft'^ by 

COS013 — — cos 91 cosflg — sin 6*^ sin03cos((p3 + ipl) . 

(A6) 

The angles are shown for illustration in Fig. S) It is 
obvious that the integration dfl as well as the one over 
either ipl or (/jg are trivial, such that Eq. (IA3|) simplifies 
to 

d$4(P;Pi,---,P4) 

-j^ I'M — ms— 7714 I" M — ^ S 1 2 

rf\/si2 / dy/s^ 



X / dcoseidcose'sdfslplWpsWql . (A7) 



Appendix B: Tensor- meson contributions to 
77, rj' 4n° 

1. Amplitude, decay width 

In this Appendix, we discuss an alternative, resonance- 
driven decay mechanism for the decays ij, r( — > 47r", 
namely, via two /2(1270) tensor mesons, see the right 
panel of Fig. [2] The Lagrangian for the tensor-tensor- 
pseudoscalar interaction reads 



TTP 



9TTP 



(Bl) 
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where the tensor nonet is given by 



\/2 \/2 



T — 



V 



K. 



(B2) 



f2^^v / 



The couphng constant gxTP is not easily determined phe- 
nomenologically; we will first write the resulting decay 
width as a function of gxTP, and then proceed to esti- 
mate it in Appendix IB 21 using QCD sum rules. 

The decay of the ^into two pseudoscalars is described 
by the Lagrangian [2l|, [s^l 

Cf, = 5t/2m.(^^"'^> = ||/2^.(9'^'^5» + . . . . (B3) 

Using the polarization sum for a tensor meson [sij 

(B4) 

with X^i, — g^i, — p^pj^/M'i , it is straightforward to de- 



rive the J2 — > TTTT decay width as [21 



r(/2 



9t 



J2 ^ 



{Mf 



AMlf\ (B5) 



Inserting Mf^ = 1275.1 MeV, T = 185.1 MeV, and 
^{f2 — >■ tttt) = 84.8% (neglecting the corresponding 
uncertainties), the coupling constant can be obtained as 
gr = 39.4 MeV. 

Applying the polarization sum of a tensor field as given 
in Eq. (jB4p . the decay amplitude for the process 77' — >■ 
72/2 47r'' takes the simple form 

/ 2 2gTTPgT V a P 

-4/2/2 = -y 3 Yi £t.^aPPlP2P3Pi 



12, S23, Sl4, S34; Si3, S24j 

- H(S13, S23, Sl4, S24; 512,534) 

- H(S12, S24, Sl3, S34; Sl4, S23)] 



(B6) 



The three terms in the square brackets are due to 
interchange of identical pions in the final state, and 



w, X, y; s, t) is given by 



'H{v,w,x,y;s,t) 



V — w — X + y 
{Ml-s){Ml^t) 



(B7) 



where we have again neglected the /2 width in the prop- 
agators. 

We calculate the partial width as in Sec. IIII C[ and 
find as the result for the decay mechanism through two 
virtual /2 states 

r(7/ ^ /2/2 ^ 4^") « 1 X 10-i6|^ MeV. (B8) 

The value of grrp is estimated using the method of QCD 
sum rules in Appendix IB 21 to be about 9GeV~^. Thus, 
the branching fraction is 



S(?7' -> 72/2 ^ 47r°) sa 4 X 10" 



14 



(B9) 



It is orders of magnitude smaller than the value in 
Eq. and hence can be safely neglected. 



2. Estimate of qttp via QCD sum rules 

In this Appendix, we estimate the unknown coupling 
constant grrp using QCD sum rules [H,!!!]. We choose 
to estimate it from the a2/27r coupling. Since we are 
not aiming at a precise calculation, complications due to 
mixing with giuon operators and anomalous dimensions 
will be neglected. The interpolating fields for the [l^l 
and TT"*" are 



j\^i^\x) = ^d{x){'j^,^u + ^u^t,)u{x) , 



= imqd{x)j5u{x) 



(BIO) 



where = (i^^ — ^)^)/2, with the standard co- 
variant derivative, and is the light quark mass. The 
fiavor wave function for the /2 is {uu + dd)/V^, and the 
corresponding interpolating field follows from the above 
equation. We will study the three-point correlation func- 
tion 



J 



n^i/a^i(p','?) 



T{j({:^)(x),f''(y)[jif'(.)]'} 



2^0 



n(p', g)p'''g'^ (.9Q^e,3i/po- + ga,^<^l3fj.pa- + gpfi<^aiypcr + gpv<^a^ip(j) + 

I 



(Bll) 



In the last step, only one Lorentz structure is kept, which The operator product expansion for the correlation 
we will use for the estimate. function can be calculated in the deep Euclidean region. 
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FIG. 5: QCD sum rule estimate of grrp- The interval 
M% G [0.8, 1.0] GeV^, marked by the vertical dashed lines, 
is optimal to suppress both higher power corrections and con- 
tinuum contribution in Ref. [34] . 



Keeping operators of the lowest dimension only, we find 



n(p',g) = 



mq{qq) 1 
16^2 



[log(V) + log(-p'')] , (B12) 



with p — q + p' , where we have neglected all polynomial 
terms since they will not contribute after the Borel trans- 
form. The correlation function can also be expressed in 
terms of hadronic quantities, which reads 



n(p',q) 



QTTpflM^ 



(B13) 

where SU(3) symmetry is assumed for the 02 and /2 
by requiring their decay constants and masses to be the 
same. The hadronic quantities are defined as 



(B14) 



1 



with 



as defined in Eq. (jB4p . Note tha t fr defined 



in this way is dimensionless. Following Ref. (ssf . we take 
p2 _ ^ perform the Borel transform only once. We 
obtain the sum rule 



Ottp 



4V^/|M6 



(B15) 



where Mb is the so-called Borel mass, and mq{qq) — 
has been used. In addition, fr was already 
calculated in QCD sum rules [s^j- Because of a cancel- 
lation between the gluon condensate and the four-quark 
condensate, the sum rule for fx is dominated by a per- 
turbative contribution, which reads ISJl 



1607r2 



-s/M 



^ds, (B16) 



where sq is a threshold parameter introduced to mimic 
the spectral function in the region > by the one 
calculated using perturbative QCD. Finally, we obtain 



grpTP 



(B17) 



This estimate is plotted in Fig. [5] as a function of il/^ . 
Taking the same interval of^M| e [0.8, 1.0] GeV^ and 
so = 2.5GeV2 as in Ref. 
estimated as 



[34| . the coupling constant is 



grprp p w 9 GeV 



(B18) 
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